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A  VEHICLE  TRAJECTORY  METHOD  FOR  INTERCEPTING 
AN  ALERTED,  EVADING  CONTACT 


1.  INTRODUCTION 


Various  military  and  nonmilitary  applications  involve  vehicle  control  systems  that  include 
functions  such  as  detection,  tracking,  classification,  localization,  vehicle  employment  (i.e.,  contact 
prediction,  vehicle  setting,  and  launch),  and  situation  assessment  (reference  1 ).  The  vehicle 
employment  portion  often  includes  the  computational  elements  for  determining  the  vehicle 
trajectory  that  will  intercept  a  designated  contact  (evader)  (reference  2),  or  place  a  vehicle 
(pursuer)  in  the  vicinity  of  that  contact,  so  that  its  internal  seeker  can  then  operate.  The  intercept 
trajectory  has  always  been  of  primary  importance  in  the  delivery  of  a  vehicle  to  a  contact,  partly 
because  this  path  results  in  the  minimum  time  to  impact.  Essentia!  parameters  associated  with  this 
vehicle  employment  problem  are 

1 .  The  deflection  angle  or  gyro  angle  (for  a  situation  where  the  pursuer  is  launched  from  a 
nontrainable  launcher), 

2.  The  time  to  initiate  seeker  search  (when  considering  vehicles  with  internal  acquisition 
and  homing  systems),  and 

3.  The  time  to  intercept  the  contact. 

The  equations  that  allow  for  the  solution  of  this  problem  form  the  basis  for  the  algorithms 
required  in  the  vehicle  control  portion  of  many  systems.  Although  the  equations  are  nonlinear, 
rapid  solution  convergence  is  required  when  used  in  real-time  applications.  Combat  control 
system  weapon-order  generation  assumes  that  the  contact  velocity  states  remain  fixed  from 
vehicle  launch  to  intercept  (reference  3).  Because  advances  in  contact  (evader)  sensor  capabilities 
have  made  it  unlikely  that  a  vehicle  would  remain  undetected  by  the  contact  until  vehicle 
intercept,  advanced  concepts  should  focus  on  the  alerted,  evading  contact  scenario. 
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Reverse  Blank 


2.  PROBLEM  DEFINITION 


The  general  problem  entails  determining  those  parameters  that  are  necessary  for 
employment  of  a  vehicle  on  an  intercept  trajectory.  For  the  acoustic  torpedo  employment 
problem,  this  requires  determining  the  straight-line  trajectory  that  results  in  the  laminar'  point  of 
the  weapon  intercepting  the  target.  Wherein  the  standard  approach  assumes  that  the  contact 
being  pursued  maintains  a  constant  course  and  speed,  the  approach  presented  in  this  report  dears 
with  a  contact  that  is  alerted  and  takes  evasive  action.  It  is  presumed  that  information,  with 
regard  to  the  expected  alertment  range  of  the  evader,  is  available,  and  that  an  alertment  strategy  is 
selected  (references  2  and  4)  A  key  aspect  of  tills  formulation  is  that  alertment  time  is  not 
required  to  be  known  a  prior  i.  The  main  torpedo  parameters  to  be  determined  include  gyro 
angle,  alertment  time,  and  intercept  time.  Other  parameters  such  as  vehicle  run  and  enable  run 
(run-to-seeker  turn-on)  are  related  and  are  readily  computed.  Provision  is  made  to  account  for 
both  evader  reaction  time  and  dynamics.  Further,  this  formulation  enables  a  new,  more 
meaningful,  error  criteria  to  be  introduced  that  results  in  the  minimum  computation  time  (number 
of  iterations)  for  particular  desired  placement  accuracy. 


3.  FORMULA  IlON/SOLUTION  OF  EQUATIONS 


FORMULATION  OF  EQUATIONS 

Figure  i  is  an  overall  block  diagram  of  the  syslein  being  considered.  Sensors  obtain 
information,  such  as  bearing,  on  the  contact  or-  target  of  interest.  This  information  is  used  by  state 
estimation  techniques  (extended  Kalman  filter,  maximum  likelihood  estimator,  etc.)  to  compute  an 
estimate  of  the  contact's  range,  course,  bearing,  and  speed.  This  contact  state  information,  along 
with  those  parameters  chosen  to  model  the  contacts  alertment  capability  and  evasion  strategy,  are 
also  provided  to  the  evading  contact  intercept  computational  unit.  This  unit  contains  the  pursuer 
evader  weapon  order  generation  (PEWOG)  models  that  determine  gyro  angle,  alertment  time,  and 
intercept  time;  and  it  is  this  unit  that  is  the  subject  of  this  report. 

A  functional  representation  of  the  evading  contact  intercept  computational  unit  is  shown  in 
figure  2.  The  basic  elements  of  the  unit  include  a  contact  model,  a  pursuer  model,  an  alertment 
model,  an  error  unit,  and  a  control  unit. 


The  contact  model  is  a  mathematical  description  of  the  contact  trajectory  from  launch  to 
intercept.  The  trajectory  is  segmented  into  two  major  components.  The  first  segment  is  from 
launch  to  the  alertment  point  (see  figure  3),  and  the  prealertment  estimates  of  contact  dynamics 
are  used  over  this  segment.  The  second  segment  is  from  alertment  to  intercept  Because  the 
contact  undergoes  a  delayed  (due  to  reaction  time)  evasion  maneuver  during  this  segment,  the 
propagation  of  its  position  to  the  intercept  point  is  based  on  both  pre-maneuver  and  post¬ 
maneuver  dynamics  (i.e.,  evasion  strategy)  as  well  as  vehicle  characteristics 
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figure  /.  Overall  Block  Diagram  of  System 


Figure  2.  Functional  Block  Diagram  of  the  Evading  Contact  Intercept  Computational  Unit 


Similarly,  the  pursuer’s  mathematical  representation  is  divided  into  the  same  two  trajectory 
components.  Because  the  pursuer  can  be  launched  from  a  platform  without  trainable  tubes  and  at 
any  deptli,  the  pursuer  has  to  transition  to  a  runout  depth  and  turn  to  the  intercept  course  Thus, 
the  propagation  of  vehicle  position  over  the  first  segment  is  a  function  of  weapon  settings  and 
characteristics  as  well  as  launcher  characteristics  During  the  second  segment,  alertment  o 
intercept,  the  pursuer  attains  optimum  search  dep  h  and  activates  its  sensor  Pursuer  position  is 
propagated  over  this  segment  based  on  weapon  settings,  characteristics,  and  sensor  parameters 


The  alertment  mode!  is  actually  a  constraint  that  is  imposed  on  both  the  contact  and  the 
pursuer  trajectories  and  is  a  functional  representation  of  the  alertment  range  (contact  detection 
capability),  in  terms  of  contact  ant!  pursuer  prealertment  parameters  Thus,  at  the  time  of 
alertment,  the  range  between  the  contact  and  pursuer  has  to  equal  the  alertment  range  for  a  valid 
intercept  solution  to  exist 

P-tTors  are  formed  between  the  computed  alertment  range  and  the  contact  detection  range, 
the  across  the  line  of- sight  positions  of  contact  and  pursuer  laminar  point  at  intercept,  and  the 
along  the  Ime-ot-sighi  positions  of  the  contact  and  pursuer  laminar  point  at  intercept  These 
errors  are  used  by  the  control  portion  to  deter  mute  the  next  set  of  updates  to  be  ted  back  for  the 
next  computational  cycle,  it  the  errors  are  not  less  ;han  the  cunveigeiu  e  uiftnia 


Figure  3.  Geometry  Depicting  Various  Points  in  Vehicle  Trajectories 


Following  the  development  in  reference  2,  the  equations  can  be  written  for  the  contact, 
pursuer,  and  aleitment  models  These  equations  take  into  account 

1  A  multiple  speed  pursuer, 

2.  The  ability  to  specify  a  seeker  tum-on  distance  or  have  one  computed  via  specification 
of  an  offset  distance  from  the  contact, 

3.  A  dive/climb  segment  to  optimum  search  depth, 

4  The  drift  of  the  pursuer, 

5  Contact  reaction  time  and  dynamics,  and 

6  The  situation  where  the  contact  is  intercepted  in  its  evasive  turn 


Definitions  of  the  variables  used  are  given  in  figures  4a  and  4b  and  the  list  of  symbols  on  page  ii 


Figure  4a.  Contact  Track 


Contact  Model 


Figure  4a  shows  the  contact  geometry  front  launch  to  intercept  Using  this  diagram,  the 
position  equations  for  the  contact  are  written  as 

sxc  :r  -Scta  sin(A)  -  Sc(tm  -  ta)  sin(A)  t  rc  cos(A)  rc  cos(A-Oc)  -  Lni  sin(A  -  0C) ,  (!) 


IYC  -Scta  cos(A)  -  Sc(tm  -  t3)  cos(A)  -  rc  sin(A)  '  rc  sin(  A  -  ()c)  -  !  ni  cos(A  -  0C)  ,  (A.) 


w  h  ere 


h ’lmc  *  in  'med  and  0C  0cm 

h  ’  -me  ■  ir.  atl^  * V  F!cdodli  ’in) 
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Figure  4b.  Pursuer  Track 


Fu.rsuer  Model 

Figure  4b  shows  the  pursuer  geometry,  including  launcher  parameters,  from  launch  to 
laminar  point  intercept  o" the  contact  Using  this  diagram,  the  equations  for  the  pursuer  can  be 
written  as 


ZXp«.P0  sm(B)  -  cos(B)  +  Rg  sin(B  -  Bg)  -  rp  cos(B  -  Bg)  rp  cos[0p  -  (B  -  Bg)] 
"  Sp(ta  •  tpc)  siiifOp  ~  (B  -  Bg)  +  (Df/2)(ta  tpC)] 

'  ^p(l-ct '  la)  s*n(0p  -  (B  -  Bg)  +  D,-(!a  -  tpC)  +  (Dr/2)(t(j  -  tg)] 

“  ^pd(‘e  "  ld)  sin[0p  (B  -  Bg)  +  ta  -  ipt-)  +  t<J  -  ta)  +  (Dr/2)(te  -  t,j)j 


S»ps(‘i  le> 


(B  •  Bg)  v  Dl(t, 


tpc* 


I>rC*d  -  *a)  1  Dr(k 


*d) 


(D,/2)(t,  -  tc)j 


(B  ■  Bg)  *  Dftta  -  tpc) 


’  D,(*d  «a>  f  I>rOe  -  »d)  *  0,'t,  •  te)J 


lYp  ~  PG  cos(B)  +  Pn  si.n(B)  +  Rg  cos(B  -  Bg)  +  rp  sin(B  -  Bg)  +  rp  sin[Gp-(B  -  Bg)] 

^p(^a  “  tpc)  cosl^p  "  (®  '  Bg  )  +  0-V2)(ta  -  tpc)l 
+  SpOd  “  *a)  ct>s[6p  -  (B  -  Bg)  +  -  tpc)  +  (Df/2)Bd  -  tg)] 

+  Spd(^e  ”  *d)  cos£®p  ■  (B  ~  Bg)  +  Df{ta  -  tpc)  +  -  tg^)  +  (Dr/2)(te  - t^)] 

+  SpS(tj  ~  te)  cost®p  ~  (B  -  Bg)  +  Dj-(ta  -  tpC)  +  Dr(ljj  -  t-j)  +  Dj-(te  -  td)  +  (Dr/2)(tj  -  te)] 

+  La  cos[0p  -  (B  -  Bg)  +  Dj^ta  -  tpC)  +  Dj^t^  -  t^l  +  Dj-(te  -  tj)  +  D^tj  -  te)]  (4) 


Alert  matt  Model 

The  alertment  model  constrains  the  range  between  the  contact  and  pursuer  to  be  equal  to 
the  contact  detection  range  at  alertment.  Using  the  prealertment  contact  and  pursuer  across  and 
along  the  line-of- sight  components,  the  constraint  equation  is  expressed  as 

Ra2  =  (P0  sin(B)  -  Pn  cos(B)  +  Rg  sin(B  -  Bg)  -  rp  cos(B  -  Bg)  +  rp  cos[Gp  -  (B  ■■  Bg)] 

"  Sp(ta  "  tpc)  sinj^Gp  -  (B  -  Bg)  +  (Dj-/2)(ta  -  tpc>]  +  Scta  sin(A)} 

+  (Rc  -  Scta  cos(A)  -  (P0  cos(B)  +  Fn  sin(B)  +  Rg  cos(B  -  Bg)  +  rp  sin(B  -  Bg) 

+  rp  sin[0p  -  (B  -  Bg)]  +  Sp(ta~tpc)  cos[0p  -  (B  Bg)  +  (Dr/2)(ta  -  tpc)]}  )  (5) 


Error  Model 

The  positional  errors  are  generated  from  the  contact  and  pursuer  across  and  along  the  line- 
of-sight  errors  at  intercept:.  Equaling  components  yields 

SXC  =  IXP)  or 

-  Scta  sin(A)  -  Sc(tm  •  ta)  sin(A)  t  rc  cos(A)  -  rc  cost  A  -  0C)  -  Lm  sin(A  0C) 

=  Pc  sin(B)  -  Pn  cos(B)  +  Kg  sin(B  -  Bg)  -  rp  cos(B  -  Bg)  +  rp  cos|Bp  -  (B  -  Bg)] 

-  Sp(ta  -  tpc)  s»n[0p  (B  Bg)  +  (D(72)(t  pc)] 

SpOd  t-a>  stnlBp  ■  (B  -  Bg)  >  (!>,/2)(td  t  ta  2t  pc)| 
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“  ^pd(^e "  ^d)  s^n[^p  “  (B  "  Bg)  +  (Df/2.)(te  +  -  2tpc)] 

'  ^ps(^i  _  *e)  s*n[®p  ~  (B  ~  Bg)  +  (Dr/2)(tj  4  te  -  2tp(_)]  -  La  sin[0p  -  (B  -Bg)  +  Df(tj  -  tpe)j,  (6) 
and  Rc  +  ZYC  «  ZYp 

Rc  -  Scta  cos(A)  -  Sc(lni  -  tg)  cos(A)  -  rc  sin{A)  +  rc  sin(A  -  0C)  -  Lm  cos(A  -  0C) 

”  P0  cos(B)  +  Pn  sin(B)  +  Rg  cos(B  -  Bg)  +  rp  sin(B  •»  Bg)  +  rp  sin[0p  -  (B  -  Bg)] 

+  SpOa  '  tpc)  cosf®p  -  (B  -  Bg)  +  (Dr/2)(ta  -  tpC)J 
+  Sp(*d  ~  *a)  cos[0p  -  (B  -  Bg)  +  (Dr/2)(td  4  ta  ~  2tpc)] 

+  Spd^e  *  *d)  cosfBp  -  (B  -  Bg)  +  (Ds/2)(te  +  t<j  ~  2tpc)] 

+  ^ps(*i  “  *e)  cos[6p  -  (B  -  Bg)  4  (Dr/2)(t;  +  te  -  2tpC)j  +  La  cos[0p  -  (B  -  Bg)  +  Dj^tj  -  tpc)].  (  7) 

Using  the  relationships  of  table  1  in  equations  (6)  and  (7),  the  positional  error  equations  are, 

L-xr  -  LxC  -  ip  cos(Bi)  +  rc  cos(A  -  0C)  +  Scta  sin(A)  +  Lm  sin(A  -  0C)  +  rp  cos(0p  -  Bi) 

-  Sp{ta  -  [Rg  4  tp(0p)]/Spj}  sin<0p  -  B,  +  (D,-/2){ta  -  [Rg  +  rp(0p).]/Spf } ) 

‘  Sp(te  -  ta  -  Ld/Spd)  sin(0p  -  Bi  +  (Dj72)(ta  4  te  -  Bd/Spd  -  2{[Rg  +  rp(0p)]/Spt})) 

-  Ld  sin(0p  -  B,  ■+■  (Dr/2)(2te  -  Lq/Spj  -  2{[Rg  4  rp(0p)j'^pt})) 

-  Sps(tj  _  ^e)  sin(0p  -  B,  +  (Dr</2)(t|  +  te  *  2{[Rg  +  rp(Qp)j/Sp£ })) 

-  La  sin(0p  -  Bj  "t  Df(tj-{[Rg  +  rp(0p)]/Sp^ } ))  ~  0  (8) 

Lyr '  Bye  4  rp  sin(B,)  -  rc  sin(A  -  0C)  f  Scta  cos(A)  +  Lm  cos(A  -  0C)  •  ip  sin(0p  -  B,) 

4  Sp{ta-[Rg  *  rp(0p)]/Spt)  cos<Op  (Bj)  *  ( f )j  / /)] ta  -  |Rg  >  ‘p(0p)i'Spj  1/ 

*  ‘’pl'e  ‘  "  Ud/Spd.)  cos(0p  -  B|  *  (1)^/2. )<,ta  t  -  I  -d/^pd  *  2 » | Rg  4  *p(0p)J/Spt })) 

+  I  d  cos(0p  -  B,  *  (lV2)<2tft  S.(|/Spp  -  2{|Rg  »  tp<()p}i/Spt  }>) 


^psv^i  ~  le)  C0Sl'^'p  '  +  (fV2){tj  +  tg  ■  2[Rg  +  rp(^p)/Sp{]  }) 

+  La  ™*<ep  -  B1  +  Dr(tj  -|R3  +  rp(©p)]/Spt>)  =  0  (9) 

In  equations  (8)  and  (9)  when  Lsto  is  selected, 

*e  "  O-sto  *  Ld)/Sp  +  Ld/Spd  , 
and  when  Lst0  is  not  selected, 

*e  ”  *i  “  WSps 

Substituting  the  relationships  of  table  1  in  equation  (5)  yields  the  alertment  range  eixor  equation, 

(Lxr  -  rp  cos(Bj)  +  rp  cos(0p  -  Bt) 

-  Sp{ta  -  [Rg  +  rp(0p)]/Sp^}  sin(0p  ••  Bi  +  (Dr/2){ta  -  [Rg  +  rp(0p)]/Spt } ) 

+  Scta  sin(A)V  +  [Rg  -  Scta  cos(A)  ••  (Lyj-  +  rp  sin(B,)  +  rp  sin(0p  -  B,) 

+  Sp{ta  -  [Rg  +  rp(Gp)]/Spt}  cos(Gp  -  B,  +  (lV2){ia  '  [Rg  +  rp(0p)]/Spt } >)] "  -  Ra2  -  0.  (10) 


1  able  I.  Sy  mbols  Defined 


Symbol 

Definition 

B, 

(B-Bg) 

b2 

(Lo/Sps  +  Ld/Spd) 

*  xc 

rc  cos(A)  -  Sc(tst)  sin(A) 

Lxr 

P0sin(B)  -  Pncos(B)  +  Rgsin(B  -  Bg) 

Lyc 

Rc  -  rc  sin(A)  -  Sc(tst)  cos(A) 

Lyr 

P0cos(B)  i  Pn  sinfB)  +  Rg  cos(B  -  Bg) 

*dive 

Ld/Spd 

le  *  *d 

(Ld)/Spd 

*m 

la  1  (st 

imc 

*-a  *  {st  1  'c^cmV^ct 

1p 

Rg/Spi 

f-pc 

[Rg  *  rp(Gp)]/Spt 

?scrh 

W-SpS 

. .  _  . . .  .  ! 

Control  Model 


The  errors  in  equations  (8),  (9),  and  (10)  are  used  to  generate  the  control  updates  required 
to  converge  to  an  intercept  solution.  The  derivation  of  the  control  is  given  in  the  next  section. 


SOLUTION  OF  EQUATIONS 

The  solution  of  equations  (8)  through  (TO)  is  now  addressed  (references  2  and  5).  The 
equations  are  transcendental  in  nature  and  do  not  lend  themselves  to  solutions  in  closed  form.  A 
numerical  solution  that  exhibits  rapid  convergence  characteristics  and  accurate  estimates  is 
required. 

Expressing  equations  (8),  (9),  and  (10)  as  general  functions  of  the  problem  unknowns  and 
performing  a  Taylor  series  expansion  yields 

eOadi.Gp)  eCtaj,tii,0pj)  +  hde/dta|j  +  j^e/a^  +  kdodGplj  +  .. .=  0  ,  (11) 

f(ta,ti,0p)  -  fltai)tii;ePi)  +  hdf/dtaij  +  jaffdtiii  +  k^aopij  + ..  ~  o ,  (12) 

g(Wi’0p)  “  8(tai>tii»0Pi)  +  h5g/5ta|{  +j5g/at;|i  +  k<3g/d0p|j  +  .  =-  0  ,  (13) 


where 


and 


ee/a.1irae/aaitana.iti=ti.i9pM)|>jj 


Ja laj  +  *i  fjj  +  j>  ®p  ’  +  k 


Neglecting  the  higher  order  terms,  the  solution  of  this  linear  set  of  equations  is  expressed  as 

h  =  [C2.2(gjC13  -  ejC33)  +  C23(e;C32  -  giC!2)  +  fj(C33C12  ~  C13C32)]/A  ,  (14) 

j  *  [fj(C3IC13  -  Cl  1C33)  +  C23(g;Cl  1  -  e,C3f)  +  C21(ejC33  -  gjCl3)I/A  ,  (15) 

k  =-  [ C22(ejC3  S  -  gjC  11)  +  fj(C32C  1 1  -  C ! 2C3 1)  -t  C2 j(gjC)2-  ejC3 2 ) )/A  ,  (16) 

and 


til 


where  t . 


Cl  l(C22C33  -  C32C23)  f  CI2(C2tC3i  -  C33C2I)  i  ('  I3«.’2K.’32  -  C3  H  22)  , 
fj,  and  g,  are  g’ven  by  equ;t*:  'ns  (8),  (9),  and  (10),  respectively 


(17) 


i  3 


The  partial  derivatives  are 


Cl  1  —  dddta  =  rc('90c/Ota)  sin(A  -  0C)  +  Sc  sin(A) 

-  cos(A  -  0C)  +  [dLm!di^  sin(A  -  0C) 

-  Sp<ta-[F.g  +  rp(0p)]/Spt>(Dr/2)  cos<0p  -  Bi  +  (D/2){ta  -  |Rg  +  rp(ep)]/Spt}> 

”  s‘nV^p  *  Bj  +  (Dp2){ta  -[Rg  +  rp(9p)j/Spt}) 

“  Sp(te  -  ta  -  Lcl/SpdXDj/2)  cos(0p  -  Bj  +  (D,-/2)(ta  +  te  -  Ld/Spd  -  2{[Rg  +  rp(0p)]/Spt})) 
+  Sp  sin(0p  -  Bi  +  (Dr/2)<ta  4-  te  -  Ld/Spd  -  2{[Rg  +  rp(Gp)]/Spt»)  .  (18) 

C12  =  dddt\  ~  rc(o0c/<9tj)  sin(A  -  0C)  -  Lm(30c/5tj)  cos(A  -  0C)  +  (<5L,n/5tj)  sin(A  -  0C) 

-  SpS(ti  -  te)(Dr/2)(l  +  dte/dti)  cos(0p  -  B(  +  (D/2X t;  +  te  -  2{[Rg  +  rp(0p)]/Spt}>) 

-  Sps(l  -  ate/5ti)  sin<0p  -  B  j  +  flV2){t}  +  te  -  2[Rg  +  rp(0p)]/Spt}> 

-  La(Dr)  cos(0p  -  Bi  +  (Dr)<tj  -  {[Rg  +  rp(0p)]/Spt}>) 

■  Sp(te  -  ta  -  Ld/Sptj)(Dr/2)(c)te/' <?tj) 

•  cos(0p  -  Bj  +  (Dr/2)(ta  +  te  -  Ld/Spd  -  2{[Rg  +  ,p(0p.)]/Spt})) 

-  Sp(5te/aq)  sin(0p  -  Bj  +  (0^2)^  +  te  -  Ld/Spd  -  2{[Rg  +  rp(0p)]/Spt}>) 

-  LdiDrXate/aj)  cos(0p  -  Bi+(DIv'2)(2te  -  Ld/Spd  2{[Rg  +  rp(0p)]/Spl}>)  .  (19) 


C 1 3  -■  ddfflp  ~  -rp  sin(0p  -  B  i ) 

-  Sp(ta  -  [Rg  *  rp(0p)]/Sp,)(  1  -  Drrp/2.Spt)  cos<0p  -  Bj  ‘  (Dr/2){ta  -  [Rg  +  rp(0p)]/Spt}> 
4  Sprp'Spf  sin(0p  -  B[  4  (D,/2 )'  tH  -  [Rg  t  *p(0p)i/Sp?}) 

'  ^p(fe  la  ”  Bd/SpdH  *  -  Dffp/Spt) 


.  cos(ep  -  B1  +  (Dr'2)(ta  +  te  -  Ld/Spd  -  2{[Rg  +  rp(0p)]/Spt}>) 

-  Ld(l  -  Drrp/Spt)  cos(9p  -  Bi  +  (Dr/2)(2te  -  Ld/Spd  -  2{[Rg  +  rp(Gp)]/Spt}>) 

"  Sps(ti  -  te)(l  -  Drrp/Spt)  cos(0p  -  B  j  +  (D/2)<ti  +  te-2{[Rg  +  rp(0p)]/Spt}>) 

-  LaO  -  IVp/Spt)  cos(0p  -  Bi  +  (Dr)<q  -  {[Rg  +  rp(0p)]/Spl}>)  .  (20) 


C21  =  dfldt&  -  rc(50c/dta)  cos(A  -  0C)  +  Sc  cos(A) 

+  LmC^c/ata)  sin(A  -  0C)  +  (dL^/dt^)  cos(A  -  0C) 

-  Sp{ta  ••  [Rg  +  rp(0p)]/Spt}(IV2)  sin(0p  -B\+  (0,7?) {ta  -  [Rg  +  rp(0p)]/Spt}> 

+  Sp  cos(0p  -  Bi  +  (D,v 2){ta  -  [Rg  +  rp(0p)]/Spt}> 

-  Sp(te  *  ta  -  Ld/Spd)(D/2)  sin(Gp  -  Bi  +  (pj2)(t&  +  te  -  Ld/Spd  -  2{[Rg  +  rp(0p)]/Spt}» 

-  Sp  cos(0p  -  Bi  +  (Dr/2)<ta  s-  tg  -  Ld/Spd  -  2{[Rg  +  rp(ep)]/Spt»)  .  (21) 


C22  -  dffdti  =  rc(39c/dtj)  cos(A  -  0C)  +  Lm(a0c/&j)  sin(A  -  0C)  +  (dLm/dtj)  cos(A  -  0C) 

‘  Sps(ti  ’  te)W2)(l  +  dtfc/dtj)  sin(Gp  -  Bj  +  (Dr/2)<tj  +  te  -  2{[Rg  +  rp(0p)]/Spt}>) 
+  Sps(l  -  dt'Jdti)  cos(0p  -B[+  (Dj72)(tj  +  te  -  2{[Rg  +  rp(0p)J/Spt ») 

-  La(Dr)  sin(0p  -  Bi  +  (Dr)<tj  -  {[Rg  +  rp(0p)]/Spt}>) 

~  Sp(te  -  ta  -  Lcj/Spcj)(Dr/2  )(7>t.c  / <9tj)  sin(0p  -  B  | 

+  (Dj72)(ta  +  te  -  Ld/Spd  -  2{[Rg  i  rp(0p)J/Spt }» 
t  Sp(cfte/t>tj)  cos(Op  -  B  i  f  (D,/2)<ta  *  tt.  L.d/Spd  2{[Rg  *  ip(0p)i/Spt}>) 

-  Ld(Dr)((^te/(>tj)  sin(0p  Bj  i  (l>r/2)<2te  Ld/ Spd  ■  2{[Rg  +•  rf)(Op>|/Spt>>)  . 


(22) 


C23  ~  df7d0p  -  tp  cos(0p  -  Bj) 


'  Sp{ta  -  [Rg  +  rp(0p)]/Spt } ( 1  -Drrp/2Spt)  siii(0p  -  Bj  +  (Dr/2){ta  -  [Rg  +  rp(0p)]/Spt}> 

■  (^prp/^pt)  cos(0p  -  B]  +  (Dr/2){ta  -  [Rg  +  *p(0p)]/Spt})-  Sp(te  -  t.a  -  L(j/Sp(j)(l  -  J-Vr^Spt) 

•  sin(0p  -  Bi  +  (Dr/2)(ta  te  -I^/Spj  -  2{[Rg  +  Jp(Qp)j/Spt})) 

-  Ld(l  ■  Drrp/Spt)  sin(0p  -  B|  t  (Dr/2)(2te  -  Ld/Sptj  -  2{[Rg  +  rp(0p)]/Sp{})) 

"  Sps(ti  -  te)(l  -  DrVSpt)  sin(0p  -  Bi  +  (D^Xq  +  te  -  2{[Rg  +  rp(Gp)]/Spt} >) 

-  La(l  -  Drrp/Spt)  sin(0p  -  Bi  +  (Dr)(tj  -  {[Rg  +  rp(0p)]^pt}))  •  (23) 

C31  -  dg/dt&  =  2{PAR1 }  [-Sp{ta  -  [Rg  +  rp(0p)]/SptKDr/2) 

•  cos(9p  -  Bi  +  (D^Xta  -  [Rg  +  rp(0p)]/Spt}> 

-  Sp  sin(0p  -  Bl  +  (Dr/2){ta  -  [Rg  +  rp(0p)]/Spt}>  +  Sc  sin(A)] 

+  2{PAR2}[-SC  cos(A)  -<-  Sp{ta  -  [Rg  +  rp(0p)]/Spt}(Dr/2) 

•  sin(9p  -  Bi  +  (lV2)(ta  -  [Rg  +  rp(0p)]/Spt}> 

+  Sp  cos<0p  -  Bi  +  (D,72){ta  -  [Rg  +  rp(0p)]/Spt }>)]  .  (24) 

C32  -  dg/dtj  -  0  .  (25) 

C33  =  dg/dQp  -  2 {PAR! }(-  rp  sin(0p  -  B \ ) 

~  SpXa  ■  [Rg  !  rp(Gp)l/Spt)(  I  -  DrJp/2Sp()  cos(0p  -  B]  t  (Dr/2){ta  -  |Rg  t  rp(0p)l/Sp^}) 
t  (Spfp/Sp,)  ssn<0p  -  Bl  t  (I-)r/2){ta  -  [Rg  l  rp(()p)j/Spt }>)  4  2{PAR2}[  -  {rp  cos(0p  -  B|) 

‘  SpUa  “  |Rg  *  rp(0p}]/,Spt }( !  ~  ?  Vi |V 2Np( )  sin(Op  -  Bj  i  ( Dr/2 )  { I  a  ■  [Rg  (  rp(0p)|/Spt}) 
(Spip/Spt)  cos<Op  -  Bj  4  (!),/2){ta  f  Kg  t-  rp{0p)|/Spt })}  j  (.26) 


where 


PARI  =  [Ljy-  -  fp  cos(B i)  +  fp  cos(0p  -  B  j) 

"  Spi^a  "  P^g  +  rp(9p)3/Spt}  sin(0p  -  Bi-KIV2){ta  ~  f^g  +  rp(^p)]^pt})  +  Scta  sin(A)] 
PAR2  =  [Rc  -  Scta  cos(A)  -  {Lyj-  +  rp  sin(B  l)  +  rp  sin(0p  -  Bj) 

+  ^p{*a "  Pg  +  fp(®p)]^pt)  cos(8p  ~  Bi  +  (Dr/2){ia  -  [Rg  +  rp(0p)]/kpt})}]  , 

and  for 

fi^mc-  ~  Sca(tj  -  Ia  -  tsj  -  Rc(^cmV^ct)*  ^nr/^a  ”  ‘^ca,  50c^a  ~  0  » 

tj  tmc.  Lfji  0,  3Ljj^/c?ta  ~  0,  and  <30g/^ta—  -0gdot , 
and  when  Lsto  is  selected 


*e  O-’sto  “  i'd)/Sp  +  Ljj/Spj,  dt^Jdi\  —  0  , 


and  when  Lsto  is  not  selected 

te  —  tj  -  Lo/SpS  ;  i^tg/Stj  —  1 


Figure  5  is  a  flowchart  of  the  PEWOG  algorithm  used  to  determine  the  estimates  for  gyro, 
alertment  time,  and  intercept  time  for  an  evading  target  scenario.  Three  logic  conditions  are 
included  in  this  block  diagram.  The  first  decision  point  addresses  the  change  in  the  computations 
dependent  on  whether  the  contact  is  intercepted  in  or  after  completion  of  the  turning  maneuver. 
The  second  decision  point  accounts  for  whether  the  seeker  turn-on  (STO)  distance  is  specified  or 
computed  based  on  a  seeker  offset  distance  relative  to  the  pursuer  position  at  intercept  The  third 
decision  point  determines  if  the  computations  have  converged  sufficiently  to  satisfy  specified 
criteria. 


Figure  5.  Flow  haft  for  Parameter  Computations 


4.  RESULTS 


Figure  6  is  a  generic  geometry  for  the  pursuer/evader  problem  in  which  the  key  problem 
parameters  are  indicated.  The  parameters  for  each  run  are  specified  in  table  2  and  were  used  by 
the  PEWOG  algorithm  to  compute  a  gyro  angle,  alertment  time,  and  intercept  time,  as  given  in 
table  3.  The  pursuer  vehicle  is  then  launched,  in  the  simulation  and  uses  the  computed  gyro  angle 
to  establish  its  course.  The  simulation  continues  to  propagate  all  vehicle  models  until  the 
specified  alertment  range  is  reached.  At  this  time,  the  contact  executes  its  evasion  maneuver, 
taking  into  account  the  reaction  time  and  contact  dynamics  specified.  All  vehicle  models  continue 
to  be  propagated  until  intercept.  The  resulting  along  and  across  the  track  errors  are  given  in  table 
3.  Figures  7  through  12  are  the  trajectory  plots  corresponding  to  the  geometries/ solutions  given 
in  tables  2  and  3. 


i’tfture  6.  iieometry  Depicting  key  Scenario  Parantetees 


1  o 


Table  2,  Run  Parameters 


Run  # 

|B 

(deg) 

A 

(deg) 

®cm 

aca/ 

^■a 

(kyd) 

ysc 

Spd^c 

^'ps^c 

STO 

Mode 

STO 

Value 

m 

-90 

90 

— 

1 

4.0 

3.1 

msm 

2.3 

Yes 

-1000 

-90 

90 

-180 

2 

4.0 

3.1 

2.7 

2.3 

Yes 

-1000 

3 

-90 

-90 

(§0 

2 

4.0 

2.3 

2.0 

2.3 

Yes 

-1000 

4 

-90 

90 

-180 

2 

■EOI 

3  1 

2.7 

2.3 

Yes 

-1000 

5 

90 

90 

45 

2 

|  4  8 

3  6 

3.3 

3.1 

No 

700 

6 

166 

Li?  . 

45 

2 

!  4.5 

3.6 

3.3 

3  1 

Yes 

-2000 

Rc  =  10,000  yards  (Runs  #1  through  #4);  Rc  =  5000  yards  (Runs  #5  and  #6),  La  =  1000  yards, 
rc  =  80  yards,  rp  =•■  128  yards,  Sct  -  Sc,  Spt  =  Sp,  I  Dr  j  =  4.3(10)"^deg/sec  (all  runs) 

Note:  Run  #1  is  a  contact  nonmaneuver  case. 


Table  3.  Solution  Parameters 


Run 


Computer  Parameters 


Gyro 


(sec) 


fci 

(yet) 


-66.7 


206.6 


324.7 


0.24 


4 

5 

6 


-73.2 


2.03  6 


99  5 


295  0 


320.3 

430.2* 


-66  9 
103.4 


215  7 


366  1 


9  5 


1  19  5 


179  3 


29,0 


1 29  0 


0.64 


-0  53 


0.86 


0.06 


i  3 


Error 


ec 

(yd) 

0  0i 


0  06 


0  01 


0  12 


0  0 1 


0  05 


Various  capabilities  of  the  pursuer  evader  tcchim|ue,  demonstrated  l>v  Ore  run  geometries 
specified  in  table  2,  are  summarized  as  follows 

1  Run  #1  treats  the  case  of  a  constant-course,  constant  speed  contact  and  is  used  as  a 
baseline  for  comparison  In  this  ran,  although  the  contact  is  alerted  at  4000  yards,  no  evasive 
maneuver  is  performed,  thereby,  resulting  in  she  solution  tor  a  nonrmmeu  voting  contact 


2.  Run  #2  uses  the  same  initial  problem  geometry  as  Run  #1;  however,  when  the  pursuing 
vehicle  closes  to  4000  yards  from  the  contact,  alertment  occurs  and  an  evasive  maneuver  is 
executed,  i.e.,  a  180°  change  in  direction  and  an  increase  in  speed.  The  gyro  angle  computed  for 
Run  #2  is  larger  in  amplitude  than  that  of  Run  #1,  reflecting  the  requirement  to  intercept  a  contact 
moving  in  the  opposite  direction  from  its  initial  course  (or  the  course  of  Run  #1). 

3.  In  Run  #3  the  contact  portion  of  the  geometry  is  a  mirror  image  of  Run  #2.  This 
demonstrates  the  ability  of  the  algorithm  to  compute  an  intercept  solution  in  a  different  quadrant 
of  the  x-y  trajectory  plane  and  for  a  contact  with  a  positive  versus  negative  turn  rate,  as  used  in 
Run  #2.  Slower  pursuer  speed  parameters  were  also  selected  for  both  the  pre-enable  and  dive 
portions  of  the  Run  #3  pursuer  trajectory. 

4.  Run  #4  shows  the  ability  of  the  model  to  compute  a  solution  when  the  point  of  intercept 
occurs  during  the  contact  evasion  turn.  The  geometry  is  the  same  as  that  of  Run  #2;  however,  the 
alertment  range  is  reduced.  This  allows  the  pursuer  to  get  closer  to  the  contact  prior  to  alertment, 
thereby,  resulting  in  the  contact  having  insufficient  time  to  complete  its  evasive  maneuver  before 
intercept  occurs.  Comparison  of  the  computed  gyro  angle  for  this  run  with  that  of  Run  #1 
(nonmaneuvering  case)  reveals  a  small  difference,  as  would  be  expected 

5.  The  PEWOG  algorithm  has  two  modes  of  operation  to  determine  the  weapon  seeker 
turn-on  (STO)  run  distance  Run  #5  shows  the  determination  of  contact  intercept  solutions  when 
a  STO  value  is  selected  as  opposed  to  being  computed  based  on  an  offset  distance  from  the 
predicted  evading  contact's  intercept  point  (as  was  the  case  for  all  other  runs)  In  addition,  for 
Run  #5,  the  initial  contact  range  at  launch,  a'ertmenl  range,  weapon  speeds,  and  launcher  tube 
(starboard)  were  selected  to  be  different  than  those  of  the  previous  runs 

6  Run  #6  examined  the  ability  of  the  algorithm  to  generate  a  solution  for  an  evading 
contact  that  requires  a  large  gyro  angle 

The  PEWOG  algorithm  takes  into  account  vehicle  drift.  Except  for  Run  i/4,  in  which  the  drift 
is  negative  (southern  hemisphere),  the  runs  are  for  a  positive  drift  rate  (northern  hemisphere). 
Examination  of  the  along  (e|)  and  across  (ec)  the  weapon  track  errors  in  table  3  reveals  the  high 
degree  of  accuracy  of  the  algorithm,  i.  e  ,  all  errors  are  less  than  1  yard  Figure  13  is  a  plot  that 
shows  the  behavior  of  the  solution  parameters  as  a  function  of  number  of  iterations  of  the 
algorithm  This  plot  is  for  Run  it?.  of  table  2,  but  it  is  characteristic  of  the  rapid  convergence 
(3  to  4  iterations)  exhibited  in  all  runs 
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Figure  7.  Trajectory  Plot  for  Run  #/ 
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higure  8  Trajectory  Plot  jar  Run  if  7 
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Figure  1 1.  Trajectory  Plot  for  Run  #5 
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figure  l  Trajectory  Plot  (or  Run  H6 
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5.  CONCLUSIONS 


The  intercept  method  developed  in  this  study  relaxes  the  current  intercept  requirement  of 
constant  contact  velocity  states  for  the  entire  intercept  run.  This  method  includes  the  capability 
for  computing  the  solution  for  a  pursuer  undergoing  multiple  speed  changes,  various  depth 
changes,  drill,  and  different  sensor  activation  criteria  In  addition,  enhancements  that  account  for 
contact  alertment,  reaction  time,  and  evasion  strategies  are  added.  Intercept  by  the  pursuer  can 
be  computed  for  any  point  in  the  contact's  trajectory,  including  during  the  portion  of  the  trajectory 
where  the  contact  is  turning  to  an  evasion  trajectory.  Finally,  the  formulation  of  the  models  and 
the  iterative  solution  technique  allow  the  intercept  parameters  to  be  computed  faster  than  real¬ 
time  (i.e.,  as  fast  as  the  previous  nonevasion  intercept  technique). 

A  high-fidelity  simulation  was  developed  to  analyze  performance.  This  computer  simulation 
includes  a  contact  vehicle  model,  launching  platform  model,  pursuer  model,  and  the  PEWOG 
equations  described  herein.  A  large  number  of  pursuer /evader  geometries  were  run,  and  the 
accuracy  and  iteration  performance  of  the  technique  was  examined.  An  illustrative  set  of  run 
geometries,  along  with  the  associated  solution  data,  are  included  to  illustrate  the  performance  of 
the  technique.  Th  e  robust  behavior  of  the  technique  is  demonstrated  by  its  ability  to  generate 
accurate  solutions  rapidly.  The  along  (ef)  and  across  (ec)  the  weapon  track  errors  between  the 
contact  and  the  pursuer  guidance  point  ,  as  determined  by  the  simulation  at  the  computed  intercept 
time,  are  also  included  along  with  a  plot  that  demonstrates  the  rapid  convergence  of  the  method 
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